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Abstract 

The resonance states of 4 H, A He and 4 Li, embedded in the three-cluster d + N + N 
continuum, are investigated within a three-cluster model. The model treats the Pauli 
principle exactly and incorporates the Faddeev components for proper description of 
the boundary conditions for the two- and three-body continua. The hyperspherical har- 
monics are used to distinguish and numerate channels of the three-cluster continuum. 
It is shown that the effective barrier, created by three-cluster configuration d + N + N, 
is strong enough to accommodate two resonance states. 
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1 Introduction 



In this paper we study the nature of resonance states in A H, A He and A Li. All these nuclei 
have a rich structure of resonance states pQ. There are 4 well-determined resonances in 4 if 
and A Li, and up to 15 resonance states were detected in A He. Most of these resonances 
have a width that is much larger than the resonance energy when measured from the lowest 
threshold. Although these resonances have been experimentally confirmed, they can hardly 
be observed in current theoretical model descriptions of these systems through standard 
elastic and inelastic scattering quantities such as S-matrix elements, differential or total 
cross sections, and so on. 

For many years, the four-nucleon system was studied by different microscopic and semi- 
microscopic methods. Different forms of the Schrodinger equation (differential ( 12112]), inte- 
gral (jU0 El) , integro-differential [7J |H1 El HH] > matrix or algebraic ( [XT] ITS} IT5] ) . ...) have been 
used to study these nuclei. Special attention was paid to 4 He, the only nuclear 4-particle 
system featuring a bound state. The theoretical study of the "ground state" of 4 if and 
4 Li, and of the excited states of all three nuclei were investigated mainly through resonance 
state analysis. Of all resonances, the first excited + state has received most attention. 
In none of the descriptions the role of the three-cluster channels was properly considered, 
and so resonance states induced by this channel could not be theoretically discovered and 
analyzed. 

Our objective is to determine the type and nature of resonance states in 4 ii, A He and 
A Li that are reproduced within a three-cluster microscopic model. For all three nuclei we 
will consider one single three-cluster configuration d + N + N. This is certainly the most 
dominant three-cluster channel, as it has the lowest energy threshold. Moreover one can 
easily construct all binary channels for these nuclei within such description. In 4 iie for 
example this configuration allows to study resonances created by the two-cluster channels 
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p+ 3 H, n+ 3 He and d + d, as well as by the three-cluster channel d + p + n. The latter should 
be very important, because 7 resonance states were detected above the d + p + n threshold. 
It is interesting to note in the same context that all four resonance states of A H lie below 
the three-cluster d + n + n threshold, whereas in A Li two resonances are found above the 
d + p + p threshold. 

We propose a modification of the method formulated in [HJ[TIiJ[Tn], an d used in [T7J[TH| to 
study resonances embedded in the three-cluster continuum, and reactions with three cluster 
exit channels. The method was shown to provide a suitable instrument for investigating 
Borromian nuclei (for instance 6 He) and nuclei with prominent three-cluster features (like 
6 Be). We wish to extend the method proposed in ^J^J to handle systems in which binary 
channels play a prominent role, by including the correct boundary conditions for both binary 
and ternary channels. The results obtained in and in [T71 ITS] allow us to restrict the 

model space to the most relevant subspace. 

To reach this objective we have to: 

• specify the microscopic modelling of the three-cluster configuration, and the approxi- 
mations to be used in calculations, 

• construct a set of dynamic equations that take into account the proper boundary 
conditions for both binary and ternary channels, 

• implement reliable numerical methods to calculate continuous spectrum wave functions 
and ^-matrix elements. 
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2 Model and Methodology 



We propose the following trial wave function for the 4-particle systems 



# = #i + #2 + #3 



(1) 



A {$i (Ax) $ 2 (A 2 ) $ 3 (A 3 ) [A (xi, qi) + f 2 (x 2 , q 2 ) + / 3 (x 3 , q 3 )]} , 



where $ a is the antisymmetric and translationally invariant internal wave function of 
the A a nucleon system, and x a , q Q are the familiar Jacobi coordinates denoting (x a ) the 
relative position of two of the clusters, and (q Q ) the relative position of the third cluster 
with respect to the former two-cluster subsystem 



with (a, (3, 7) a cyclic permutations of (1, 2, 3). The three components of the wave functions 
{^1, ^2, ^3} (more precisely {fi, f 2 , fz}) have to be determined by solving the many-particle 
Schrodinger equation. Specific symmetries of the system can reduce the number of compo- 
nents: if the three-cluster configuration contains two identical clusters, only two distinguish- 
able components {fi, f 2 } will occur; this is the case for 4 if and 4 Li. If all three clusters 
are identical (impossible for the 4-particle system though), only one independent component 
{fi} would occur. 

We shall use a matrix or algebraic form of the Schrodinger equation. To this aim we 
expand the wave function f a (x a , q a ) in an oscillator basis (referred to as a BiOscillator 
(BO) basis): 





a ' 



(2) 
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where 

\n y ,l,n x ,X;LM) a = ® nyjl (q a ) $ nxjX (x a ) {Yi(q a ) ■Yx(jc a )} LM (3) 
and <3> ni ; (q) is the familiar (radial) oscillator function: 

The total angular momentum L of the three s-clusters is a vector sum of two partial angular 
momenta l a and X a associated with the respective Jacobi coordinates x Q and q a . 

As each cluster function $ a (A a ) is antisymmetric by construction, the antisymmetriza- 
tion operator in (JTJ) only involves the permutation of nucleons between clusters, and it can 
be represented as 

A = 1 + A12 + A23 + A 31 + A 123 , (5) 

where A a p exchanges nucleons from clusters a and (3, and A123 permutes particles from all 
three clusters. In some respects this antisymmetrization operator is similar to a short range 
interaction. Indeed, the antisymmetrization is noticeable only when the distance between 
clusters is small. At larger distances both the potential energy and the antisymmetrization 
effects become negligibly small. The operator .4.123 is important only when the distance 
between all three clusters is less than a certain restricted value. If one of the clusters (say, 
cluster a) is far away from the two other clusters (3 and 7, and the latter are close to one 
other, Apy will have a pronounced contribution, as well as the two-cluster interaction Vg 7 
derived from the iViV-potential. 

Each set of Jacobi coordinate x a and q a , and each set of oscillator functions (j3J) 

{\n y ,l a ,n x ,X a ;LM) a } 

for a = 1, 2 or 3 cover the whole configuration space (i.e. account for all possible relative posi- 
tions of three clusters in space). We will limit ourselves to the subspace {\n y , l a = L, n x , X a = 0; LM)} 
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of the total Hilbert space. Two arguments for such a choice can be given (in particular for 
the 4-particle system): 

1. We deal with s-shell clusters, and the two-cluster compound subsystems also are s-shell 
nuclei; the latter (d, 3 H and 3 He) have a ground state containing a dominant S-wave 
contribution. 

2. It was shown in P^ITH IT^] that this subspace dominates in the full Hilbert space. For 
instance, the ground state energy of 6 He and 8 He obtained within this subspace differs 
less than 1% from the energy obtained in the full Hilbert space. It was also shown that 
this subspace dominates in the wave function of the 2 + state of e He, appearing as a 
resonance embedded in the three-cluster continuum a + n + n, as well. 

To include the proper boundary conditions, we will split the oscillator space ^C^J lrix A | 
(a = 1, 2, 3) into the internal and asymptotic parts. The former consists of the basis functions 
of the lowest iVj oscillator shells (i.e. all functions with N = n y + n x = 0, 1,2,... ; 
it involves (iVj + 1) (iVj + 2) /2 basis functions for each value of a). With this size of the 
internal region in oscillator space one can evaluate the maximal size of the three-cluster 
system in coordinate space by using the correspondence between oscillator and coordinate 
representations (see details in [201121]) 

iCax ^ by/4Ni + 6. (6) 

If, for instance, the oscillator length b = 1.5 fm and Ni = 15, then -R max — 12 fm. 
The maximal distance between any pair of clusters will be of the same magnitude. -R max 
also corresponds to the minimal size of the three-cluster system in the asymptotic region. 
So in the internal region all three clusters are close to each other, which means that all 
antisymmetrization components are important, as well as all interactions between clusters. 

6 



In the asymptotic region we distinguish two different regimes. In the first regime the 
distance between two clusters is small, while the third one is far apart. In the second regime 
all three clusters are well separated. If a selected pair of clusters (say, j3 and 7 clusters) has 
(a) bound state (s), then the first regime is responsible for scattering of the third cluster on 
the compound j3 + 7 subsystem. This process can be described by two-body asymptotics. 
The second regime is associated with the full disintegration of the three-cluster system, 
with three independent (non-interacting) clusters. These two regimes have to be treated 
differently. This means that two different forms of the wave function have to be used to 
properly describe these processes. It will require some reconstruction of the basis functions 
in order to suit both two- and three-body physical processes in the exit channels. 

In the first regime of the asymptotics we can neglect all antisymmetrization components 
but A/3 7 . As for the potential energy, the most important contribution is generated by the 
two-cluster potential V^ 7 . The other components V a p and V aj , originating from the short- 
range iViV-forces, are negligibly small, and only long-range Coulomb forces are of importance. 
The wave function in coordinate and oscillator representation can then be factorized as 



The two-cluster wave function and its counterpart in oscillator space < C„ A >, is an 



* a = A {$1 (A) $2 (A 2 ) $3 (A 3 ) f a (x a , q Q )} 

~ Afr (Ap) % (A^) g a (x Q )} $ Q (A a ) / W (q a ) 

= *i 2) • $ a (A a ) / W (q Q ) 

= E C ^i,n x ,x\ny,l,n x ,\;LM) a 
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eigenfunction of the two-cluster Hamiltonian H a 

h^= E f * + E 

By solving the Schrodinger equation 

A'2 



(7) 



E { n *> x 



H (2) _ E (*) 







n x =0 



with a chosen number of basis functions N2, we obtain the bound state(s) (a = 0, 1, . . .) 
of the two-cluster subsystem, which determine the threshold energy of the two-body break 
up of the tree-cluster system. 

In the second regime of the asymptotics we can neglect the antisymmetrization operator 
and the short-range components of the inter-cluster potential. In this regime, we use the 
Hyperspherical Harmonics (HH) basis to describe the full decay of the three-cluster system, 
because (see, for instance, [221 1231 EH Ej ) this basis is the obvious choice for such type of 
three cluster behavior. The transition from the bioscillator basis \n y , I, n x , A; LM) a to the HH 
basis \n p , K; I, X; LM) a (see details of the definition of HH functions in e.g. [15]) is performed 
by an orthogonal matrix. This transformation can be calculated in a straightforward way. 

The asymptotic part of the wave function will then be represented by two sets of expansion 
coefficients 



n y ,L 1 °n p ,ii" min ;L' ( - / n p ,ii" min +2;L' 



where K m i n = L. All expansion coefficients in the asymptotic region have a similar form, and 
consist of incoming (v4 + \ Wk ) anc ^ outgoing (ijr L ~\ ip^ ) waves (see detail of the definition 
in H3IIBI): 



u n s ,L — 



2-Rr, 



3c ;a ^ (j^aRny) S Co -, a lfj\^ {k a R ny ) 



(+) 



C n p ,K;L ~ ( 2 PnY K^K ^ K ' {kpn) ~ S Co , aK ( k Pn 



H 



,(+) 



(10) 

(11) 
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where the index cq denotes the entrance channel and 



2m 



k=\j —E, p n = b^jAn p + 2K + Q 



k 



2m 



E — E. 



(«) 



. - \i -p- v j - J -n R n y = b^jAn y + 2L + 3. 



(12) 
(13) 



Note that the index a numerates the binary channels, while both indexes a and K distinguish 
the ternary channels. Thus Co equals ao, if the entrance channel is a binary one, or Co = «o> Kq 
for the three-cluster entrance channel. 

With this definition of the asymptotic part of the wave function we deduce the equations 
for the scattering parameters and wave function. By taking into account (JBJ), Q, (JTUJ) and 
(ITU) one obtains 



a n y ,n x <Ni 



5 n p >Ni 
5 n y >Ni 
5 n„>Ni 



H -E 

H — E n 
H-E 



n y ,n x )__ C~J~ x 



a \ flyi Tlx 



4 5 t ^ (*si2B.) 



^3/ > 



a \ i ^x 



H-E 
H-E 



y,E^\ ^ {k~ a R~ ny ) 
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pi 



(14) 



This system of linear equations contains in the general case of three different clusters 



- (AT- + 1) (Ni + 2) + 3 + 3N ch . HH 



(15) 



parameters to be determined. Here 



Nch.HH = (K max — K min ) /2 + 1 



(16) 



is the number of HH channels. From this total amount, | (iVj + 1) (iVj + 2) coefficients 
represent the wave function in the internal region, and the other 3 + 3N c h.HH parameters 
determine the elastic and inelastic processes, leading to two or three clusters in the exit 
channels. These parameters unambiguously define the wave functions in the asymptotic 
region. 

3 Results 

We use the Minnesota (MP) j23j, and the modified Hasegawa-Nagata (MHNP) [23 EE] 
nucleon-nucleon (NN) potentials. The oscillator radius b is chosen to optimize the bound 
state energy of the deuteron, and equals b = 1.512 fm for MP and b = 1.668 fm for MHNP. 
Considering these two potentials reveals the effect of peculiarities of iViV-forces on the pa- 
rameters of resonance states. 

In a first calculation, we neglect all binary channels, and only consider the three-cluster 
channels. This allows to understand what kind of resonances are generated by this channel 
only. We have omitted spin-orbital components of the iViV-forces to reduce the compu- 
tational burden. Results obtained in this approximation can be considered to represent a 
"lower limit" for the width of a resonance, as additional channels will open new decay pos- 
sibilities of the resonance, which will increase its width. In this respect the three-cluster 
approximation will indicate whether some resonances state (s) could survive after binary 
channels are included. 

By solving the dynamic equations ()14j) for N c channels, we directly obtain the N c x iV c 
^-matrix. There are two different parametrizations for the S'-matrix. In the first one, 
each element can be represented by the phase shift 8^ and the inelastic parameter T]ij: 
Sij = rjij exp {2i5ij}. In the second parametrization the S'-matrix is reduced to diagonal form 
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by an orthogonal transformation: 



S 



O 



T 



S (E) 



o 



(17) 



The latter procedure leads to N c uncoupled elastic "eigenchannels" , whose (eigen)phase shifts 
parametrize the diagonalized S^'-matrix. We use the eigenphase shifts to determine both 
the energy and width of the resonances. They are obtained from the following equations 



We start our analysis from eigenphase shift. In Fig. Q we display eigenphase shift of 
so-called 3=^3 scattering for L 71 = l~, 5=1 state in the 4 Lz, obtained with the Minnesota 
potential. Similar pictures are obtained for other nuclei and different (L n , S) states and also 
with the MHN potential. One can see from Fig. Q that there are two resonance states in 
4 Li, first resonance is narrow and manifest itself through the first eigenchannel, while the 
second resonance is very broad and appear in the third eigenchannel. 

The eigenphase shifts provide with a direct and simple way to determine the energy and 
width of a resonance state in the three-cluster continuum, but to get information concerning 
the main features of the three-cluster dynamics one needs to analyze the phase shifts 5ij and 
the inelastic parameters rjij. 

In Figs. |2Jand|21we display the phase shifts 5u and inelastic parameters rja connected with 
the diagonal matrix elements of the original S-matrix. They are obtained for A Li and with 
MP. One notices that only two hyperspherical harmonics are responsible for the lowest l - 
resonance state. The K = 3 phase shift displays the classical behavior for a resonance state, 
while the K = 5 phase shift indicates a "shadow" resonance. Many more hyperspherical 
momenta are involved in creating the second resonance. 

In Tables ^and El we collect the parameters of the resonance states lying above the three- 
cluster threshold d + N + N. The even parity states are obtained with i^ max = 10, and the 
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Figure 1: Eigenphase shifts for the L w = 1 , S — 1 state of 4 Li obtained with the Minnesota 
potential and K max = 11. The eigenchannels are numerated. 
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Figure 3: Inelastic parameters of the 3=^3 scattering for the L w = 1~, S = 1 state in 4 Li. 
Results are obtained with the Minnesota potential and K max = 11. 



14 



odd parity states with K ma , x = 11. 



Nucleus 


L w 


S 


K 

1 v max 


E, MeV 


T, MeV 


E, MeV 


r, MeV 


4 H 


\~ 





11 


1.642 


0.367 


6.726 


2.759 


4 H 


\- 


1 


11 


1.911 


0.374 


6.958 


2.982 


4 Li 


\- 





11 


2.604 


0.413 


7.787 


3.141 


4 Li 


\- 


1 


11 


2.912 


0.465 


8.085 


3.384 


4 He 


2+ 





10 


1.950 


0.233 


2.904 


0.207 



Table 1: Resonance states of 4 if , A He and 4 Li, created by the three-cluster channel d+N+N. 
Results are obtained with the Minnesota potential. 



Nucleus 


U 


S 


K 


E, MeV 


T, MeV 


E, MeV 


T, MeV 


4 H 


\- 





11 


3.972 


1.170 


9.469 


3.440 


*H 


\- 


1 


11 


3.738 


0.950 


9.250 


3.362 


4 Li 


\- 





11 


0.748 


0.093 


5.009 


1.531 


4 Li 


\- 


1 


11 


0.662 


0.056 


4.772 


1.329 


4 He 


2+ 





10 


0.890 


0.005 


2.436 


0.167 



Table 2: Resonance states of 4 if, 4 He and 4 Li, created by the three-cluster channel d+N+N. 
Results are obtained with the MHN potential. 

It is known that there is an effective barrier in each channel of three-cluster system. The 
barrier is created by a potential well, resulted from a NN-interaction between nucleons from 
different clusters, and centrifugal barrier, which is proportional to 

h 2 K(K + 4) 
2m p 2 
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In 4 He and A Li the Coulomb repulsion 

Z eff 
P 

have to be added to the effective barrier. The effective charge Z e ff depends on quantum 
numbers a, K, li, I2 and its definition can be found in [IS] . The deeper is potential well, 
the larger is the effective barrier and, consequently, more resonance states can be created 
by such effective barrier. One can see that the effective barrier, generated by the MHN 
potential, is more deeper than the one connected with the Minnesota potential. As a result 
of this difference, the resonance states in A H, obtained with the MHN potential, have larger 
energy then resonances, calculated with the Minnesota potential. 

As we pointed out the results, presented in Tables [T] and |21 for the even parity states 
are obtained with K max = 10, and for the odd parity states with K max = 11. These values 
for K max are sufficient to obtain stable results. This is demonstrated in Fig. H] where the 
parameters (energy and width) of the 1~ resonance in 4 H are displayed as a function of 
if max . The results in Fig. 0] are presented for the Minnesota potential, and similar results 
are obtained for the MHN potential. We indicate some "false" resonances states appearing 
at small values of K max due to the restriction on decay channels compatible with this K max . 
When we increase the number of open channels, the "false" resonances disappear and the 
physical resonances converge to their correct positions. 

There are some arguments that the physical resonances do not depends on the boundary 
conditions implemented, or on the used approximations. One can e.g. study the behavior 
of the so-called Harris states, i.e. the eigenvalues of the Hamiltonian as a function of the 
number of basis functions involved in a calculation. It was shown (see, for instance, [27| I28|) 
that the eigenvalues E v (n) (E v (n) is u-th eigenvalue, obtained with n basis functions) 
create plateaus at the energies of resonance states. For a very narrow resonance this plateau 
is already observed for a small value of n. For wider resonances, one needs more basis 
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max 



Figure 4: Energy of the 1~ resonance state of 4 if (total spin S — 1) as a function of K max . 
Error bars indicate the double of the resonance width. The calculations have been performed 
with the Minnesota potential. 
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functions to reach a plateau. For a small number of basis functions a wider resonance can 
become apparent as the repulsion of two eigenvalues (avoided crossing of two eigenvalues). 
Such behavior of these eigenvalues was observed for the Hamiltonian of the three-cluster 
configuration d + N + N in 4 H, 4 He and 4 Li. 

4 Conclusion 

A microscopic model is formulated to treat properly the two- and three-body boundary 
conditions. For this aim the Faddeev component is used. The hyperspherical harmonics are 
used to numerate three-cluster channels. They are very valuable for describing three-cluster 
asymptotic. Two AW-potentials are involved in the calculations in order to evaluate the 
sensitivity of the final results with respect to this important factor of a microscopic model. 

The model is applied for studying resonances states in 4 H, 4 He and 4 Li nuclei, created 
by the three-cluster configuration d + N + N. The results presented here demonstrate 
that the three cluster configuration creates an effective centrifugal barrier which allows to 
accommodate several resonances. The effect of the two-cluster channels on the position and 
width of the three-cluster resonances will be discussed in future work. 
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